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Abstract: We reveal that reahstic fluids generate microscopic- level discontinu- 
ity constantly and the discontinuity spreads out with motion of particles rather 
rapidly and widely. These things cannot be treated by the standard kinetic 
equations, and thus the existing fluid theories, macroscopic ones and micro- 
scopic ones, need to be revised considerably. 



It has been unanimously believed that the following standard framework | 
Continuity equations and Liouville's theorem 
i 

BBGKY hierarchy equations 

. ^. . 

Kinetic equations — > Fluid equations 

i i 
Equilibrium and nonequilibrium phenomena 
gives sound microscopic and macroscopic descriptions of all classical fluids ob- 
served in the physical reality. 

Nevertheless, discoveries in various fields, such as those related to turbulence, 
chaos and dissipative structures, constantly revealed fluid phenomena that ap- 
pear to be inconsistent with our existing theoretical knowledge. In the studies 
of statistical behavior of charged particlesQQ, we ourselves were in trouble to 
construct intuition out of what was obtained from the kinetic equations. The 
situation motivated us to make a lot of investigation on related fundamental 
subjects and finally resulted in our critical attitude toward the existing kinetic 
theory. 

In the usual textbook elaboration of the kinetic theory[^0], all the particles 
are at first assumed to make motion in the way as if they constitute a contin- 
uous medium and the mediumlike behavior is cast into a differential operator 
similar to that in fluid mechanics; an integral operator, as a correction term, 
is then attached, in which particles interact with each other as if they are or- 
dinary colliding particles. In our view, the formalism, though seems to give 
a comprehensive description of the two contradictory aspects, medium aspect 
and particle aspect, runs into conceptual difficulties. According to the theory, 
if the number of particles becomes more and more (the gas becomes denser 
and denser), then the coUisional term becomes overwhelming; otherwise, the 
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medium aspect becomes more important. (In some textbook, it is explicitly 
stated that if the collisions between particles are not significant the Boltzmann 
equation is reduced to the coUisionless Boltzmann equation, called sometimes 
the Vlasov equation.) However, this picture is not in harmony with the simple 
physical intuition which says that a "fluid" consisting of a very small number 
of particles behaves as a collection of particles only. The particles in such a 
"fluid" move like individual particles and collide with each other like individual 
particles; no mediumlike character can be observed. 

In this letter, we will investigate the issue thoroughly. It will be shown 
that a realistic Boltzmann gas (by "Boltzmann gas" we mean a gas in which 
interforces between particles are short-ranged) is capable of generating several 
kinds of discontinuous distribution functions and the dynamical behavior of such 
distribution functions is far from mcdiumlike and cannot be described by the 
existing kinetic theory. 

We start our discussion by recalling the standard kinetic equations briefly. 
These equations can be written in the following unifying form 

collision 



dt dr m dv J \ St 



where the right side represents an integral coUisional operator, which takes 
different forms for different kinetic equations. For purposes of this paper, we 
confine ourselves to a gas in which collisions are not significant. Under this 
understanding, Eq. (|l|) becomes 

d d F d \ ^ ^ 

ot or m ov J 

Moving together with one particle in the fluid, we simply find that the equation 
can be expressed as a convective derivative 



dt) 



= 0. (3) 

path 



In words, this is to say that the distribution function keeps invariant along any 
particle path in the six-dimensional phase space (/i-space) . 

The invariance expressed by (^ is quite amazing in the sense that the dis- 
tribution function can have such nice and regular behavior while each particle 
in the fluid seems free to make any kinds of motions, including some "irregular 
ones" . This is by no means possible and seeking for hidden requirements is 
kind of necessary. A careful inspection leads us to the following three require- 
ments, (i) Throughout the dynamical process, the distribution function must 
be continuous both in terms of position and in terms of velocity, (ii) The force 
term F in the equation must be free from dissipation, (iii) No stochastic forces 
(collision-type) get involved, namely the force F must be sufficiently smooth in 
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space. In what follows, we are mainly concerned with the first requirement and 
it will be shown that a realistic fluid can constantly generate discontinuity at 
the miroscopic level and the discontinuity spread out rather rapidly and widely. 

Figure 1 




(a) (b) 

Firstly, consider a piece of boundary shown in Fig. 1. On the right side 
of it, there are particles belonging to an ordinary gas initially. Suppose that 
we are able to observe these particles continuously. It is easy to find that the 
distribution function at every point on the left side of the boundary involves 
discontinuity in terms of velocity angle. 

Figure 2 




Then, let's look at the particles leaking through a small hole of a gas con- 
tainer (free-expansion gas), shown in Fig. 2. We can find that the distribution 
function related to all these particles is continuous in the position space but not 
in the velocity space. 

Figure 3 




Actually, there are more sources of discontinuity. To get a complete knowl- 
edge about this, it is instructive to study another kind of interaction between 
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particles and boundaries, collisions between particles and boundaries. Such col- 
lisions are usually assumed to be kind of deterministic and the assumption is 
adopted by many college- level textbooks , in which molecules are considered as 
having hard-spheres, walls as perfectly smooth and collisions as classically elas- 
tic, shown in Fig. 3. According to this elastic model, the distribution function 
must keep invariant after the collisions, which partly explains why boundary 
problems have not received much attention in the standard kinetic theory. 

Figure 4 
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It is easy to see that the model outlined above is an ideal one serving for 
educational purposes. Actual observations in fluid mechanics tell us a different 
story [Q: the fluid particles immediately next to a solid surface remain almost 
stationary with respect to the surface (apparently in the average sense), shown 
in Fig. 4. This is called the no-slip condition by scientists in fluid mechanics. 

Figure 5 




To conform to the no-slip condition, a different model, called the statistical 
model, should be adopted. In this model, particles leave, after collisions, the 
solid surface in the way as if they are emitted from there, and their velocities 
will distribute according to a statistical law irrespective of how they come to 
there. (Though more sophisticated models may be better off in the academic 
sense, we take this simple one as our "zeroth-order approximation" for relative 
simplicity.) 

Taking this statistical model implies that we are ready not only to admit 
that interaction between particles and boundaries is of a dissipative nature and 
stochastic nature but also to admit that the distribution functions of emitted 
particles are discontinuous: if the area of the surface is small enough the situ- 
ation is like the gas leaking out of a container in Fig. 2; if the surface is of a 
finite-size, the distribution function will be roughly the same as that created by 
the finite boundary in Fig. 1. 
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Examples aforementioned have shown that reahstic boundaries can indeed 
generate, by blocking particles or by emitting particles, certain kinds of dis- 
continuous distribution functions. A conception which may come to one's mind 
immediately is that by introducing the (5-function and/or the step function these 
discontinuous distribution functions may be treatable as if they are continuous 
ones. Such conceptions got a lot of successes in other physical fields, why not in 
this field? We will see, however, that the position- velocity phase space, in which 
distribution functions are defined, possesses special features and one of the fea- 
tures is that the position r and the velocity v are assumed to be completely 
independent of each other. (Note that we have v = r for each particle in the 
fluid.) It is because of this assumption that the existing kinetic theory is inap- 
plicable in the situations even if the 5-function and step-function arc properly 
introduced. 

We investigate the following distribution function 

/=^'5(v-z;o^), (4) 

where uq and Vg represent two constants, r the position vector and r = |r|. By 
referring to Fig. 2, we can recognize that the distribution function characterizes 
a kind of free expansion gas. The time development of these particles can be 
expressed by 

x' =^ X + Vxt 

y' = y + vyt (5) 

z' ~ Z + Vzt 

and 

Vx=Vx, v'y=Vy, v'z=V,_. (6) 

By a direct calculation, we know that along a particle trajectory 

/(^,r,v) ^ 



fit' 

Actually, Fig. 2 can directly tell us that, in contrast with (||), the distribution 
function is not invariant along the particle trajectories. 

The difficulty becomes even more striking if we substitute the distribution 
function ^ into the kinetic equation (|^) . The following two partial derivatives 

d[6{v-vor/r)] 



dr 

and 

=? (9) 



d[6{v — vor/r)] 



i9v 

are rather puzzling. It is obvious that the troublesome situation arises not 
because we do not know how to differentiate the ^-function but because we 
cannot determine whether or not r and v correlate. 
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Up to now, we have seen that such discontinuous distribution functions pose 
difficulty in the two related respects: (i) The discontinuity spreads throughout 
almost the entire space instead of "staying" on certain boundaries, (ii) Since the 
discontinuity spreads the distribution function must contain path information 
of particles; this is in conflict with the partial-differential equations where r 
and V are considered independent of each other. The latter respect, called the 
paradox of r-v correlation in this letter, holds its significance rather generally. 

Before finishing this letter, it is important and relevant to look at whether 
or not other related theories, such as the continuity equations and Liouville's 
theorem, hold in the presence of such discontinuity. (In fact, by referring to 
the paradox of r-v correlation, readers themselve can come up with the answer 
to the question without many studies.) Adopting that the force term in (H) is 
independent of velocity, we find that the coUisionless Boltzmann equation is 
identical to the continuity equation in /j,-space 

df d(vf) d(vf) „ 
ot or ov 

It is then obvious that the continuity equation in /i-space does not generally 
hold for discontinuous distribution functions. 

In a similar way, it can be shown that the continuity equation in the grand 
phase space (F-space), which can be regarded as a natural extension of the 
continuity equation in /i-space, is not valid for systems having discontinuous 
distribution functions. 

Finally, we study a more fundamental issue: Is Liouville's theorem applicable 
in statistical mechanics? To get a definite answer to the question and, at the 
same time, to gain insight into the important subject, it is necessary to examine 
derivations of the theorem in detail. 

In textbooks, one of the derivations is based on the continuity equation in 
F-space and will be discussed no more. Another derivation|^, which invokes 
the Hamiltonian formalism and the Jacobian approach, is worth our special 
attention. 

Suppose that the Hamiltonian of a system with 2N degrees of freedom is in 
the form 

^ = E^ + *(9i'92. •••.*)• (11) 

N 

We have explicitly assumed that the potential term $ only contains the position 
coordinates <Zi , 92 , • • • and implicitly assumed that the potential term $ is differ- 
entiable in terms of these coordinates. Note that these two things have defined 
the syatem as a pure mechanical one with no dissipative and stochastic forces 
involved. Suppose that dT represents a phase volume element in F-space. One 
version of Liouville's theorem states that the value of dT is motion-invariant. 

To prove it, we adopt that the time-development of the system is a canonical 
transformation 

p{to),q{to)'^Pit),Q{t). (12) 
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In the grand space we have, generaUy, 



dqdp = DdQdP, (13) 

where 

5(0, P) 



is the Jacobian of the transformation. By applying the properties of the canon- 
ical transformation, we get (with the details omitted) 

D = l. (15) 

Equation ( |l3| ) can then be written in the form 

dToito) = dVit). (16) 

Eq. (^6|) is regarded as Liouville's theorem in classical mechanics. To have 
Liouville's theorem in statistical mechanics, which states that the grand density 
p conserves in statistical processes, several requirements have to get involved. 
By no coincidence, the requirements are almost exactly the same as those we 
listed for the invariance of distribution functions: In the entire statistical process 
the grand density p has to keep continuous and the external forces, which are 
associated with $ in (^ij), must be completely free from the dissipative nature 
and stochastic nature. As has been shown, these requirements cannot be gen- 
erally fulfilled for realistic statistical systems. All these tell us that Liouville's 
theorem should be regarded as a theory in classical mechanics only. 

In summary, this letter has revealed that seemingly simple fluids may actu- 
ally be beyond the scope of the existing fluid mechanics and beyond the scope 
of the existing kinetic theory. 
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